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Abstract. In this paper we study how to approximate the Leray weak 
solutions of the incompressible Navier Stokes equation. In particular 
we describe an hyperbolic version of the so called artificial compress- 
ibility method investigated by J. L. Lions and Temam. By exploiting the 
wave equation structure of the pressure of the approximating system we 
achieve the convergence of the approximating sequences by means of dis- 
persive estimate of Strichartz type. We prove that the projection of the 
approximating velocity fields on the divergence free vectors is relatively 
compact and converges to a Leray weak solution of the incompressible 
Navier Stokes equation. 



This paper is concerned with the convergence of the artificial compress- 
ibility approximation to the Leray weak solutions ("turbulent in the Leray 
terminology" ) of the 3-D Navier Stokes equation on the whole space. This 
approximation was introduced by Chorin [HO], Temam |3UU31| and Oskolkov 
|21| , in order to deal with the difficulty induced by the incompressibility con- 
straints in the numerical approximations to the Navier Stokes equation. The 
paper of Temam |3U1 13 lj and his book [32] discuss the convergence of these 
approximations on bounded domains by using the classical Sobolev com- 
pactness embedding and they recover compactness in time by the classical 
Lions ^21 method of fractional derivatives. This paper will take a differ- 
ent point of view, namely we wish to exploit the underlying wave equation 
structure and the presence of dispersive type estimates. In particular we 
will consider the following system 



where (x,t) G R 3 x [0,T], u £ = u £ (x,t) G R 3 and p e = p £ (x,t) G R, 
p = f 6 (x,t) G R 3 . 



The system will be discussed as a semilinear wave type equation for the 
pressure function and the dispersive estimates will be carried out by using 
the L p -type estimates due to Strichartz . The particular type of 

Strichartz estimates that we are going to use here can be found in the book 
of Sogge [57j or deduced by the so called bilinear estimates of Klainerman 
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1. Introduction 



d t u £ + Vp £ = /j,Au e 
edtp 6 + div u £ = 0, 




(1.1) 
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and Machedon |14j and Foschi Klainerman jHJ. Our analysis can also be 
related to the convergence of the incompressible limit problem via a formal 
expansion (see for instance Temam |32) . Chapter 3). In particular a simi- 
lar wave equation structure has been exploited in various way by the paper 
of P.L. Lions and Masmoudi ^Bl, Desjardin, Grenier, Lions, Masmoudi [I], 
Desjardin Grenier [JJ]. 

In this paper we analyze the convergence problem in the case of the whole 
space but our method can be extended to exteriors domains which will be 
done in a forthcoming paper. 

The interest into the artificial compressibility methods started with the pre- 
viously mentioned results of Chorin and Temam and was later on investi- 
gated by Ghidaglia and Temam |2j . Later developments of numerical inves- 
tigations in the directions of projections methods have been carried out by 

PH, H> E2, 120], E2], M- 

This paper is organized as follows. In Section 2 we recall the mathematical 
tools needed in the paper and recall same basic definitions. In Section 3 
we set up our problem, we explain our approximating system and we state 
our main result. The Section 4 is devoted to recover the a priori estimates 
needed to get the strong convergence of the approximating sequences and to 
prove the main theorem. Finally in Section 5 we give the proof of the main 
result. 

2. Preliminaries 

For convenience of the reader we establish some notations and recall some 
basic facts that will be useful in the sequel. 

We will denote by V(R d x R + ) the space of test function Cg°(R d X R+), 
by V'(R d x R+) the space of Schwartz distributions and (•, •) the duality 
bracket between V and V and by M t X' the space C°([0, T];X)'. Moreover 

W k ' p (R d ) = (I - A)~^L p (R d ) and H k (R d ) = W k > 2 (R d ) denote the non- 
homogeneous Sobolev spaces for any 1 < p < oo and k G R. W k,p (R d ) = 
(I-Ay^L p (R d ) and H k (R d ) = W k ' 2 (R d ) denote the homogeneous Sobolev 
spaces. The notations L P L% and L P Wx ,q will abbreviate respectively the 
spaces L p ([0,T];L q (R d )), and L p ([0, T];W k,q (R d )). 

We shall denote by Q and P respectively the Leray's projectors Q on the 
space of gradients vector fields and P on the space of divergence - free vector 
fields. Namely 

Q = X7A- 1 div P = I-Q. (2.1) 

Let us remark that Q and P can be expressed in terms of Riesz multipliers, 
therefore they are bounded linear operators on every W k)P (1 < p < oo) 
space (see |28])- 

Let us recall that if w is a (weak) solution of the following wave equation in 
the space [0, T] x R d 

(-g + A^w(t,x)=F(t,x) 
w(0,-) = f, d t w(0,-)=g, 
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for some data f,g,F and time < T < oo, then w satifies the following 
Strichartz estimates, (see "TO], [15] ) 

\H\l\li + ll^llif^ £ WSWm + hh-r 1 + I^ILfrg" ( 2 - 2 ) 

where (q,r), (q,r) are wave admissible pairs, namely they satisfy 

and moreover the following conditions holds 

1 d d Id 

- + - = tj-7 = - + --2. 

q r 2 q' r 

Later on we shall use (|2.2|) in the case of d = 3, (q',r') = (1,3/2), then 
7 = 1/2 and (q,r) = (4,4), namely the following estimate 

Ml*. + \\dtw\\ Lfw -iA < \\f\\ A x/2 + 11511^-1/2 + ll*1l £ i £ 3/a. (2-3) 

Beside the Strichartz estimate (|2.2|) or ()2.3|) in the case of d = 3 (see [2*7]). 
a more refined estimate, related to an earlier linear Strichartz [22] estimate, 
can also be deduced by the bilinear estimates of Klainerman and Machedon 
[14] . Foschi and Klainerman [""], namely 

Ml*. + \\d t w\\ Lfw -x,4 < \\f\\ A i /2 + \\g\\^ /2 + \\F\\ LjL 2. (2.4) 

3. Approximating system and main result 
Let us consider the incompressible Navier Stokes equation 

dtu + div(u (8) u) — fiAu = Vp + / 
divu = (3.1) 
u(x,0) = n , 

where (x,t) £ M 3 x [0, T], u £ R 3 denotes the velocity vector field ,p£R 
the pressure of the fluid , / £ M 3 is a given external force, fi is the kinematic 
viscosity. Let us recall (see P. L. Lions [19] and Temam "5""]) the notion of 
Leray weak solution. 

Definition 3.1. We say that u is a Leray weak solution of the Navier Stokes 
equation if it satisfies (|3.1|) in the sense of distributions, namely 

n( Vn ■ Vip — UiUjdiipj — u ■ -^- \ dxdt 



(f, ^H-^mdxdt + / n • <pdx, 

for all tp e Cf?(R d x [0, T]), divip = and 

divu = in V'(R d x [0,T]) 
and the following energy inequality holds 



1 



\u(x, t)\ 2 dx + \i I I | Vn(x, t)\ 2 dxds 
JO JR d 



< — f \uo\ 2 dx+ [ (f,u) H -i xH ids, for alH > 0. 
2 jRd J 
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There exists in the mathematical literature several results concerning the 
existence of Leray weak solutions to the Navier Stokes equations, for example 
we can refer to books of P.L. Lions |19| and Temam [32]. The case d = 3 
is a major open problem and a considerably more difficult case than the 
case d = 2, since the bound on the L 2 norm (kinetic energy) provides only 
a control on a supercritical norm and does not provide any information 
concerning the critical controlling (and scaling invariant) norm L 3 . Hence 
we do not know (opposite to the case d = 2) whether or not the Leray weak 
solutions are unique, unless (see Serrin [21]) we assume a control on the L 3 
norm. Some important regularity results can be found in pp. 
In order to approximate the system 1)3. ljl we wish to use the system 1)1.1)1 
where we introduce a "linearized" compressibility constraint given by the 
equation 

1 

dtp £ = — divu £ . 
e 

In order to avoid the paradox of increasing the kinetic energy along the 
motion we introduce the correction 

1 

— -(divu £ )u £ 
into the momentum balance equation. 

The limiting behaviour as e J of the initial data to 1)1.1)1 deserves a little 
discussion. Indeed (jl.l)) requires two initial conditions 

u £ (x, 0) = u £ {x), p £ {x, 0) = p £ (x), (3.2) 

while the Navier Stokes equations require only one initial condition on the 
velocity u. Hence our approximation will be consistent if the initial datum 
on the pressure will be eliminated by an "initial layer" phenomenon. Since 
in the limit we have to deal with Leray solutions it is reasonable to require 
the finite energy constraint to be satisfied by the approximating sequences 
(u £ ,p £ ). So we can deduce a natural behaviour to be imposed on the initial 
data (uq,Pq), namely 

ul = n e (-,0) — ► n = u(-,0) strongly in L 2 (R S ) (ID) 

Ji v l = ve/(-,0) — ► strongly in L 2 (R 3 ). 

Let us remark that the convergence of \fsp% to is necessary to avoid the 
presence of concentrations of energy in the limit and it includes the Temam's 
assumption that {pq} is bounded in 1? . 

Since it will not affect our approximation process, for semplicity from now 
on, we will take fj, = 1 and f £ = 0. For convenience, let us now formulate 
an existence theorem concerning the approximating problem ()1.1)1 . 

Theorem 3.2. Let (uq,Pq) satisfy the conditions (ID) for some e > 0. Then 
the system 1)1.1)1 has a weak solution (u £ ,p £ ) with the following properties 

(i) u £ € L oo ([0,T];L 2 (lR 3 ))nL 2 ([0,T];ij 1 (M 3 )). 
(ii) ^p £ ei M ([0,T];I 2 (R 3 )), 

for all T > 0. 

The proof of this theorem will be omitted since it will be a consequence of 
all the "a priori bounds" that will be obtained in the sequel and it will follow 
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from the use of standard finite dimensional Galerkin type approximations. 
Let us now state our main result. The convergence of {u £ } will be described 
by analyzing the convergence of the associated Hodge decomposition. 

Theorem 3.3. Let (u £ ,p £ ) be a sequence of weak solution in M 3 of the 
system assume that the initial data satisfy (ID). Then 

(i) There exists n G L°°([0, T]; L 2 (R 3 )) n L 2 ([0, T]; i? 1 ^ 3 )) such that 

u £ ^-u weakly in L 2 ([0, T]; i^QR 3 )). 

(ii) The gradient component Qu £ of the vector field u £ satisfies 
Qu £ — > strongly in L 2 ([0, T]; L p (M 3 )) , for any p G [4,6). 

(iii) The divergence free component Pu e of the vector field u £ satisfies 

Pu £ — >Pu = u strongly in L 2 ([0, T]; L 2 oc (R 3 )). 

(iv) The sequence {p 6 } will converge in the sense of distribution (more 
precisely in H^W* 2 ' 4 + M t Wx 1,4/3 + L^H' 1 ) to 

p = A~ x div ((u ■ V)«) = A~Hr((Du) 2 ). 

Moreover u = Pu is a Leray weak solution to the incompressible Navier 
Stokes equation 

P(d t u -Au+(u- V)u) = in V'([0, T] x M 3 ), 

and the following energy inequality holds 

\( \u(x,t)\ 2 dx+ [ [ \Vu{x,t)\ 2 dxdt<\ [ \u(x, 0)\ 2 dx. (3.3) 

Remark 3.4. This theorem can be easily extended to the nonhomegeneous 
equation (|3.1j) . by assuming 

f £ — / strongly in L 2 ([0, T]; iJ-^R 3 )). 

Remark 3.5. Let us denote by Rj the Riesz transform. The Hardy space 
W 1 (M 3 ) is a closed subspace of L X (1R 3 ) defined by 

H 1 ^ 3 ) ={/ € L 1 ^ 3 ) | Rjf G L 1 (M 3 ),for any j = 1....3}. 

Then one has 

p G Z^flO, T]; L 3 / 2 (M 3 )) n ^([0, T]; L 3 (R 3 )), (3.4) 
and there exits c\ > 0, such that 

\\{ tr ( Du ) 2 )\\Lmo,T];HHR 3 )) - c ^\\ u o\\h(w 3 y 

4. A PRIORI ESTIMATES 

In this section we wish to establish the priori estimates, independent on 
e, for the solutions of the system which are necessary to prove the 

Theorem 13.31 We will achieve this goal in two steps. First of all we will 
recover the a priori estimates that come from the classical energy estimates 
related to the system (jl.ip . Then we get stronger estimates by exploiting 
the structure of the system. In fact, as we will see later on , the sequence p £ 
satifies a wave type equation. This will allow us to apply to p £ the Strichartz 
estimates (|'2,4j) . (|2.2|) . and to get in this way dispersive bounds on p £ . 
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4.1. Energy estimates. The next results concerns the energy type esti- 
mate for the system (|1,1|) . 

Theorem 4.1. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (jl.lj) . Assume that the hypotheses (ID) hold, then one has 

E(t)+ [ [ \Vu £ (x,s)\ 2 dxds = £7(0), (4.1) 
Jo Jr 3 



where we set 



LV)- I [\\u*{x,t)\ 2 + £ -\p £ {x,t)\*)d,. (1.2, 



Proof. We multiply, as usual, the first equation of the system (|1.1|) by u £ 
and the second by p £ , then we sum up and integrate by parts in space and 
time, hence we get (|4.1|) . □ 

Corollary 4.2. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (|1.1|) . Let us assume that the hypotheses (ID) hold, then it 
follows 

y/sp £ is bounded in L°°([0, T]; L 2 (M 3 )) ; (4.3) 

ep\ is relatively compact in i/ _1 ([0, T] x M 3 ), (4.4) 

Vu e is bounded in L 2 ([0,T] x M 3 ), (4.5) 

u £ is bounded in L°°([0, T]; L 2 (M 3 )) n L 2 ([0, T]; L 6 (R 3 )), (4.6) 

(u £ -V)u £ is bounded in L 2 ([0, T]; L 1 (M 3 )) n L x ([0, T]; L 3 / 2 (R 3 )), (4.7) 
(divu s )u e is towndec/ mL 2 ([0,T];L 1 (M 3 ))nL 1 ([0,T];L 3/2 (IR 3 )). (4.8) 

Proo/. (|43|) . ijPj) . (|4~5|) follow from (gUJ), while (|4~6|) follows from (jlTTjl 
and Sobolev embeddings theorems. Finally (|4.7|) and (j4.8j) come from (|4.5|) . 

(gSJ). □ 

4.2. Pressure wave equation. In this section by using the Strichartz es- 
timates ()2.3|) . (12. 4j) we get a priori estimates on p e . We will use a wave 
equation structure for p £ . First of all let us rescale the time variable, the 
velocity and the pressure in the following way 

r = — u(x,t) = u £ (x,^t), p(x,t) = p £ (x,^t). (4.9) 

As a consequence of this scaling the system (jl.lj) becomes 

d T u + \feVp = \/eAu — \fe {u-V)u— ^j- (divu)u ^ ^ 

\fed r p + e?zt> u = 

then, by differentiating with respect to time the equation (j4. lOjt ^. and by 
using (14.10^ . we get that p satisfies the following wave equation 

d TT p — Ap + Adivu — div ((u ■ V) u + ^(divu)u J =0. (4-11) 
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Now we consider p = pi + p2 where pi and p2 solve the following wave 
equations: 

J d TT pi — Api = —Adivu = F\ 
|pi(x,0) = d T pi(x,Q) = 0, 

d T Tp2 — Ap2 = div ( (u ■ V) u + -{divu)u 



(4.12) 



2 V 7 / (4.13) 
jfa(x, 0) = p(x, 0) d T p2(x, 0) = d T p(x, 0). 

Therefore we are able to prove the following theorem. 

Theorem 4.3. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (|l.ljl . Assume that the hypotheses (ID) hold. Then we set 
the following estimate 

£3/ Vll4w t - 2 ' 4 + £7/8 W d tP £ \\ L iw t - 3 > 4 ~ V^Mhl + \\divu £ \\ H -i + VT\\divu £ \\ 

+ || {u £ ■ V)u £ + \{divu £ )u £ \\ TlT 3/2. 

(4.14) 



Proof. Since pi and p2 are solutions of the wave equations (|4.12j) . (|4,13|) . 
we can apply the Strichartz estimates ()2.3j) and (|2.4j) . with (x,r) G R 3 x 
(0, Tj^fe). Since A _1 pi satisfies the equation 

MA" 1 !*) " A(A- X p!) = A- 1 ^, (4.15) 

then by using the Strichartz estimates (|2.4j) we get 

IIA-^H^ + \\d T A~ x pi\\ L%w -iA < HA-^IUi^ (4.16) 

namely 



T 

Pillow" 2 - 4 + 11^1^4^-3,4 < j^WdivuW^. (4.17) 
In the same way we have that A~ l / 2 p2 satisfies the equation 

d TT (A-y%) - A(A~ 1 / 2 p 1 ) = A- l ' 2 F 2 , (4.18) 
therefore by using the estimate 1)2. 3|) we get 

\\&- 1/2 P2\\l* <x + \\d T A-Wp2\\L*W-V Z \\^- 1,2 p{xM\ H ^ 

+ \\A- l l 2 d T p(xM\ H -^ 

+ ||A- 1 / 2 F 2 |L 1 ,8/ji (4.19) 

namely 

H^IU^-M + \\d T p 2 \\ L ^ w -2A < \\p{x,G)\\ H -x,2 + \\d T p(x, 0)|| H - 8 /a 

+ || (u ■ V) u + -(diuu)€t|| rl r 3/ 2 (4.20) 
2, l t l x , 
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Now by using (|4.17|) . Q4.2UJI it follows that ft verifies 

\\P\\ L 4 W - 2 ' 4 + \\ d rPh± w -3A < \\fti\\ Liw -2A + \\ft 2 \\ Liw -i A (4.21) 

+ \\drftl\\ r 4 W -3A + \\d T ft2\\ rAlv -2A 



< 



\\ft(x,0)\\ H - 1/2 + \\d T ft(x,0)\\ H - 3/ 



v T 1 

+ jifi \\divu\\ L 2 L 2 xi + || (fi • V) u + - (dwfiJfiH^s/a. 

Finally, since 

II ~ll — l/2r II £ || 

\\P\\L^{{0,T/^);L<i{M?)) = £ \\P \\L r ([0,T\;L^(R 3 )) 

we end up with ()4.14|) . □ 

5. Strong convergence 

In this section we conlcude the proof of the Theorem IM.M1 In particular 
we will show that the gradient part of the velocity Qu £ converges strongly 
to 0, while the incompressible component of the velocity field Pu £ converges 
strongly to Pu = u, where u is the limit profile as e [ of u E . 

5.1. Strong convergence of Qu £ and Pu 6 . We start this section with 
some easy consequences of the a priori estimates established in the previous 
section. 

Corollary 5.1. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (|l.lj) . Assume that the hypotheses (ID) hold. Then, as e J 0, 
one has 

ep £ — >0 strongly in L°°([0, T]; L 2 (R 3 )) n L 4 ([0, T]; W^ 2 ' 4 (IR 3 )) ; 

(5.1) 

divu £ — ► strongly in W -1,oo ([0, T\; L 2 (R 3 )) n L 4 ([0, T]; W- 3A (R 3 )). 

(5.2) 

Proof. (jSHJ), (jS"2)) follow from the estimates (|Q)) . flUS)) and the second 
equation of the system (|l.lj) . □ 

Now, we wish to show that the gradient part of the velocity field Qu £ goes 
strongly to as e j 0. As we will see in the next proposition, this will be a 
consequence of the estimate (|4.14j) and of the following auxiliary result. 

Lemma 5.2. Let us consider a smoothing kernel ifj G Cfi°(M. d ), such that 
"0^0; JjRd = 1, and define 

Then for any f £ if 1 (R d ), one has 

Wf-f* ^ahpm < Cp« 1-(7 llv/L2 (Md) , (5.3) 

where 

pe[2,oo) ifd = 2, p€[2,6] if d = 3 and a = d(---). 
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Moreover the following Young type inequality hold 

11/ * ^allLp(Rd) < Ca~ d («~p) WfWw-s,^, (5.4) 
for any p,q G [1, oo], q < p, s > 0, a G (0, 1). 

Proposition 5.3. Lei us consider the solution (u e ,p £ ) of the Cauchy prob- 
lem for the system p.ljl . Assume that the hypotheses (ID) hold. Then as 
e | 0, 

Qu e — ► strongly in L 2 ([0,T}; L P (R 3 )) for any p G [4,6) . (5.5) 



Proof. In order to prove the Proposition 15. HI we split <5m £ as follows 

where is the smoothing kernel defined in Lemma 15.21 Now we estimate 
separately J\ and J2. For J\ by using (|5.3f) we get 



1 1 



Jx<a 13 U p; (/ \\VQu £ (t)\\l2dtj < a ^ \\Vu £ \\ L 2 L 2. (5.6) 

Hence from the identity Qu £ = — E 1 > 8 'VA~ 1 £ 7 ' 8 dtp and by the inequality 
(|5.4|) we get J2 satisfies the following estimate 



J^e^nvA-V/s^^n < e V8 a -2-3U-|j|| £ T/8 atp 11 



< e^a-'-^-^T^lle^dtpll^-sA. (5.7) 

Therefore, summing up Q5.6JI and (|5.7|) and by using (|4.5|) and 1)4. 14j) . we 

conclude for any p G [4, 6) that 

\\Qu £ \\ L 2 L p < Ca~ 3 (?~p) + Ct£ 1/8 cT 2 ~ 3 (^) . (5.8) 

Finally we choose a in terms of e in order that the two terms in the right 
hand side of the previous inequality have the same order, namely 

a = E 1 ' 1 \ (5.9) 

Therefore we obtain 

\\Qu £ \\ L 2 L v < C T e^ for any p G [4,6). 

□ 

It remains to prove the strong compactness of the incompressible compo- 
nent of the velocity field. To achieve this goal we need to recall here, the 
following theorem (see |25j). 

Theorem 5.4. Let be T C L p ([0, T];B), 1 < p < oo, B a Banach space. T 
is relatively compact in L p ([0, T];B) for 1 < p < 00, or in C([0,T]; B) for 
p = 00 if and only if 

(i) I J f(t)dt, f G b\ is relatively compact in B, < t\ < ti < T, 

(") I™ \\f(x + h)- f(x)\\ L pQo,T-h};B) = uniformly for any f G T. 
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The compactness can be obtained by looking at some time regularity 
properties of Pu £ and by using the Theorem 15.41 but before we need to 
prove the following lemma. 

Lemma 5.5. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (jl.lj) . Assume that the hypotheses (ID) hold. Then for all 
h € (0, 1), we have 

\\Pu £ (t + h)- J Pu £ (t)|| L 2 ([0iT]xR a ) < C T h l >\ (5.10) 

Proof. Let us set z £ = u £ (t + h) — u £ {t), we have 

\\P U £ (t + h)- Pu £ (t)\\ 2 L2{%T]xR3) = J J dtdx(P Z £ ) ■ (PZ £ - PZ £ * Va) 

+ / / dtdx{Pz £ ) ■ (Pz £ * V> Q ) =h + h- 
JO 

(5.11) 

By using (|5.3j) we can estimate I\ in the following way 

h < \\Pz £ \\ LTLl £ \\Pz £ {t) - (Pz £ * ^ a )(t)\\ L 2dt 

< aT^llul^jVul^. (5.12) 
Let us reformulate Pz £ in integral form by using the equation 1)1.1^ . hence 

rT i- rt+h 

h < dt dx ds(Au £ - (u £ ■ V)u £ - -u £ (divu £ )(s,x) ■ (Pz £ * ip a ){t,a 
Jo Jr 3 Jt 2 

(5.13) 



Then integrating by parts and by using Q5.4j> . with p = oo and q = 2, we 
deduce 



h < h\\Vu £ \\ 2 2 + Ca- z/2 T l l 2 \\u £ \\ LaaL 2 (h [ t+ \ (u £ ■ V) u £ - -{divu £ )u £ \\ 2 L 
t,x t * \ J t 2 

< h\\Vu £ \\ 2 L2 + Ca-^T^hW^W L?Ll \\ (u £ -V)u £ - l -{dtvu £ )u £ \\ L2Ll . 

t,X t X £ t x 

(5.14) 

Summing up I\, I2 and by taking into account (|4.5|) . (j4.6j) . (|4,7j) . (|4.8|) . we 

have 

\\Pu £ (t + h)- ^ £ (i)||i2 ([0jT]xR 3 ) < CiaT 1 ' 2 + ha^T 1 ' 2 + h), (5.15) 
by choosing a = h 2 / 5 , we end up with (|5.1U|) . □ 

Corollary 5.6. Let us consider the solution (u £ ,p £ ) of the Cauchy problem 
for the system (|l.lj) . Assume that the hypotheses (ID) hold. Then as e \ 

Pu £ — ► Pu, strongly in L 2 (0, T; L 2 oc (R 3 )). (5.16) 

Proof. By using the Lemma 15.51 and the Theorem 15.41 and the Proposition 
Owe get (t5~16|) . □ 
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5.2. Proof of the Theorem 13.31 (i) It follows from the estimate (|4.6[) . 

(ii) It is a consequence of the Proposition 15.31 

(iii) By taking into account the decomposition u £ = Pu £ + Qu £ , by the 
Corollary 15.61 and the Proposition 15.31 we have that 

Pu e — ► u strongly in L 2 ([0, T]; Lf oc (R 3 )). 



(iv) Let us apply the Leray projector Q to the equation then it 

follows 

Vp £ = AQu £ - Q (div(u £ ® u £ ) + -u £ divQu £ ) . (5.17) 



2 

Now by choosing a test function <p £ H}W% A/Z n C^Wl A n LfH^ and by 
taking into account (|4.5|) . (j5.5j) . (|5.1fij) . we get, as e { 0, 

(u £ divQu £ ,Q v ) < \\Qu £ \\ L 2 Lt \\Vu £ \\ L 2 L 2\\Qv\\ LrL 4 

+ \\Qu £ \\L2Lih £ \\L~L2\\VQv\\ L ?Li -» 0, (5.18) 

(div(u £ <g> u £ ),Qip) = (div(Pu £ <g> Pu £ ), Qip) + (div(Qu £ <g> Qu £ ), Q(p) 

+ (div(Pu £ Qu £ ),Qip) + (div(Qu £ <g> Qu £ ),Q^p) 

-> (dro(Pu (8) Pu),Q<p) = (Qdiv((Pu ■ V)Pu),<p). 

(5.19) 

So as e | we have, 

{Vp £ , <p) — ► (VA -1 dii;((ii • V)ti), if). (5.20) 

Now we can pass into the limit inside the system (jl.ljl and we get u satisfies 
the following equation in V'([0,T] x M 3 ) 

P(d t u - An + (u ■ V)«) = 0. (5.21) 

Finally we prove the energy inequality. By using the weak lower semicontinu- 
ity of the weak limits, the hypotheses (ID) and denoting by x the weak-limit 
of \fep £ , we have 

f ^-\x\ 2 dx+ [ \\u(x,t)\ 2 dx + / / \Vu(x,t)\ 2 dxdt 

JR3 2 J K 3 2 J J R 3 

<liminf(/ i|n £ (x,t)| 2 dx+ / % £ \ 2 + [ [ \Vu £ (x,t)\ 2 dxdt 

\Jr3 2 JR3 2 Jo iM 3 

= liminf / - (|u^| 2 - e|p§| 2 ) da; = / \\u \ 2 dx. (5.22) 
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